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ABSTRACT 

We search for correlations between the spin in pairs of spiral galaxies, to study if 
the angular momentum gain for each galaxy was the result of tidal torques imprint 
by the same tidal field. To perform our study we made use of a sample of galaxy 
pairs identified using the Sloan Digital Sky Survey. We find a weak, but statistically 
significant correlation between the spin magnitude of neighbouring galaxies, but no 
clear alignment between their orientation. We show that events such as interactions 
with close neighbours play an important role in the value of the spin for the final 
configuration, as we find these interactions tend to reduce the value of the A spin 
parameter of late-type galaxies considerably, with dependence on the morphology of 
the neighbour. This implies that the original tidal field for each pair could have been 
similar, but the redistribution of angular momentum at later stages of evolution is 
important. 

Key words: galaxies: formation - galaxies: fundamental parameters - galaxies: in- 
teractions - galaxies: statistics - galaxies: structure - cosmology: observations. 



1 INTRODUCTION 

A fundamental issue in the study of galaxy formation and 
evolution is the internal angular momentum, from its acqui- 
sition to identifying the role it plays in determining the final 
structure of galactic systems. 

To explain the origin of the internal angular momentum 
in galaxies, Hoyle (1949) proposed the tidal field theory, later 
studied by Peebles (1969) and White (1984). In this theory, 
galaxies acquire their angular momentum by gravitational 
torquing exerted, in the early Universe, by neighbour pro- 
togalaxies. The tidal fields from the large-scale structure 
surrounding the protogalaxy exert a torquing moment in- 
creasing the angular momentum of the system, which grows 
linearly in time prior to the gravitational collapse of the 
structure as protogalaxies have large linear sizes, and thus, 
the surrounding material is able to imprint angular momen- 
tum efficiently. At later stages, after the protogalaxy decou- 
ples from the expanding background and turns around, the 
growth of angular momentum is reduced to a second order 
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effect, because the collapse dramatically reduces the lever 
arms (Schaefer 2009). 

Several analytical and numerical studies have been per- 
formed to study this theory. The first concern was to quan- 
tify the amount of angular momentum imprinted on the 
galaxies (Peebles 1969; Doroshkevich 1970; White 1984; 
Barnes & Efstathiou 1987). It is usual to express the to- 
tal angular momentum of the galaxies through the dimen- 
sionless angular momentum parameter A, which in standard 
cosmological models peaks at values around Ao ~ 0.04 and 
shows a log-normal distribution with a spread in the loga- 
rithm of 0.48 < o"a < 0.66 (for a compilation of theoretical 
estimates on Ao and a\ see Shaw et al. 2006). In consistence 
with what was found in the indirect studies of Syer, Mao 
& Mo (1999), the estimates of the distribution of A for cos- 
mological dark matter haloes from analytical calculations 
and N-body simulations, have recently been shown to be 
in agreement with inferences of this quantity from observed 
galactic structural parameters for large volume limited sam- 
ples, by Hernandez et al. (2007) and Berta et al. (2008). 

Since the introduction of the A parameter, and with the 
boom of high resolution N-body simulations and detailed 
semi- analytical models, a well stocked collection of studies 
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has emerged, some looking at the evolution of the spin pa- 
rameter (Warren et al. 1992; Bullock et al. 2001; Vitvit- 
ska et al. 2002; Pierani, Mohayaee & de Freitas Pacheco 
2004, Davis & Natarajan 2009), dependencies with environ- 
ment (Avila-Reese et al. 2005; Maccio et al. 2007; Bett et 
al. 2007) and correlations with virial mass (Barnes & Efs- 
tathiou 1987; Cole & Lacey 1996; Tonini et al. 2006; Bett 
et al. 2007, Cervantes-Sodi et al. 2008), and others concern- 
ing the direct influence of this parameter on the internal 
structure of the galaxies (Fall & Efstathiou 1980, Mo, Mao 
& White 1998, van den Bosch 1998, Jimenez et al. 1998, 
Prantzos & Boissicr 2000, Boissier et al. 2001, Hernandez 
& Cervantes-Sodi 2006, Berta et al. 2008, Cervantes-Sodi & 
Hernandez 2009). Most of this works deal only with numeri- 
cal simulations, lacking observational counterparts. This, as 
the angular momentum of the dominant dark haloes can not 
be directly measured in real galaxies, although some recent 
efforts have begun to address the issue of obtaining large ob- 
servational samples as counterparts (e.g. Berta et al. 2008; 
Cervantes-Sodi et al. 2008). The orientation of the angular 
momentum, however, is easier to obtain, if one assumes that 
the shape of the galaxies is largely coupled to the overall an- 
gular momenta of the haloes. 

The interest concerning possible alignments between 
galaxies due to the mechanism of acquisition of angular mo- 
mentum, has motivated several observational studies. Gott 
& Thuan (1978), exploring the angular momentum in the 
local group, pointed out that if a pair of galaxies (like the 
Galaxy and M31) is formed in isolation, and the collapse 
phase of formation was short compared with other dynam- 
ical time scales, the spin vectors of the galaxies should be 
perpendicular to the separation vector between the galax- 
ies and parallel between them. This, for the particular case 
of the local group, holds rather well. Sharp, Lin & White 
(1979) extended the study to a sample of nearly 100 pairs 
and reported a complete lack of correlation between the 
spins. Oosterloo (1993) reached the same conclusion, con- 
trary to the result obtainned by Helou (1984), who found 
an asymmetric correlation; they report that the spin vec- 
tors avoid being parallel, in favour of the anti-parallel con- 
figuration. More recently, Pestana & Cabrera (2004) studied 
the relative alignment of spin axes for observed pairs of spi- 
ral galaxies, and found some significance against the null 
hypothesis of random orientations of binary spiral galaxies. 
They conclude that complex and repeated interactions prob- 
ably occur in binary spiral galaxies, so that spins that are 
parallel, antiparallel or nearly orthogonal often occur. Slosar 
et al. (2009) detected a correlation in the spin directions of 
pairs of spiral galaxies, with significant correlation at small 
separations (< 0.5Mpc/h). 

With the aim of elucidating the process of galaxy as- 
sembly, some studies concerning alignments of galaxies in 
clusters have also been performed. It is believed that galax- 
ies and clusters are formed through a process of hierarchical 
accretion. In that case, primordial alignments could quickly 
be erased by dynamical interactions (Coutts 1996; Plionis 
et al. 2003), and the current alignments would be produced 
by recent accretions episodes (Faltenbacher et al. 2005). The 
positive detection of alignments seems to depend on several 
characteristics of the clusters, such as the presence of sub- 
structure (Plionis & Basilakos 2002), its accretion history 
(Faltenbacher et al. 2005), the morphology of the galaxies 



(Faltenbacher et al. 2007; Torlina, De Propris & West 2007; 
Wang et al. 2009) and the morphology of the cluster itself 
(Aryal & Saurer 2006; Aryal, Paudel & Saurer 2007), which 
reflects not only the initial conditions, but the dynamical 
interactions taking place in such complex systems. 

Using numerical simulations, Barnes & Efstathiou 
(1987) fixed their attention on this issue, and found no con- 
sistent indication of coherence in the spins of adjacent ob- 
jects. The result of Bailin & Steinmetz (2005) is similar, 
they could not find a clear tendency for the principal axes 
of neighbouring haloes to point in a preferred direction, and 
Porciani, Dekel & Hoffman (2002) argue that spatial corre- 
lation of spins on scales larger than a few Mpc, induced by 
primordial tidal torques, are strongly affected by non-linear 
effects. 

In this paper we extend the study of angular momen- 
tum correlations between neighbouring galaxies, using the A 
parameter to account for the magnitude of the spin of ob- 
served galaxies. This adds an extra variable to the analysis 
of angular momentum in observed galaxy pairs, which had 
up to now been limited to considering only the orientation of 
the galaxies. A more quantitative description of the problem 
is hence possible, allowing for a more detailed comparison 
to the outcome of cosmological simulations, and yielding in- 
teresting constraints on the origin of spiral galaxy pairs. As 
a first test we begin with the simplest case, that of galaxy 
pairs. This is a regime which makes our work complementary 
to existing studies in clusters e.g. Faltenbacher et al. (2005) 
or Wang et al. (2009). Extensive studies on the effects of the 
large-scale environment and interactions with neighbouring 
galaxies on various other galaxy properties are presented by 
Park, Gott & Choi (2008, hereafter PGC) and Park & Choi 
(2009). 

The plan of this paper is as follows. In Section 2, 
we review the derivation of the A spin parameter for in- 
ferred haloes of spiral galaxies as developed in Hernandez 
& Cervantes-Sodi (2006), in Section 3 we present the SDSS 
sample of galaxy pairs used in the present work. Section 4 
presents the results from the analysis, and finally, in Section 
5, we summarize our general conclusions. 



2 ESTIMATION OF THE SPIN FROM 
OBSERVABLE PARAMETERS 

To quantify the magnitude of the angular momentum, we 
employ the A spin parameter as introduced by Peebles 
(1969); 



where E, M and L are the total energy, mass and angu- 
lar momentum of the configuration, respectively. In Hernan- 
dez & Cervantes-Sodi (2006) we derived a simple estimate 
of total A for dark halos hosting disc galaxies in terms of ob- 
servational parameters, and showed some clear correlations 
between this parameter and structural parameters, such as 
the disc to bulge ratio, the scale height and the colour. Here 
we recall briefly the main ingredients of this model. The 
model considers only two components for galaxies, a disc for 
the baryonic component with an exponential surface mass 
density E(r); 
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E(r) = S e" 



r/Rd 



(2) 



where r is a radial coordinate and Eo and Rd are two 
constants which are allowed to vary from galaxy to galaxy, 
and a dark matter halo having an isothermal density profile 
p(r), responsible for establishing a rigorously flat rotation 
curve Vd throughout the entire galaxy; 



p(r) = 



— f-V 

4ttG V r J 



(3) 



We assume: (1) that originally the baryonic and dark 
matter components of the protogalaxy were well mixed, 
aquiring the same specific angular momentum, as often 
asumed in galactic formation models (e.g. Fall & Efstathiou 
(1980); Mo et al. (1998)), and recently confirmed by sev- 
eral authors addressing this particular issue (Zavala et al. 
(2008); Croft et al. (2008); Piontek & Steinmetz (2009)); 
(2) the total energy is dominated by that of the halo which 
is a virialized gravitational structure; (3) the disc mass is 
a constant fraction of the halo mass F = Md/Mn- These 
assumptions allow us to express A as 



A = 



2 1 ' 2 V d 2 Rd 
GMh 



(4) 



Finally, we introduce a baryonic disc Tully-Fisher (TF) 
(e.g Gurovich et al. 2004) relation: M d = CWVd ' 5 , and tak- 
ing the Milky Way as a representative example, we evaluate 
F and Ctf to obtain 



A = 21.8 



-Rd /kpc 



(Vd/kms- 1 ) 3 /' 2 ' 



(5) 



In Cervantes- Sodi et al. (2008) we tested the accuracy 
of our estimate of A, comparing the value obtained using 
equation [5] to the actual value of A from numerically sim- 
ulated galaxies from six distinct groups, where the actual 
value of A is known and where we had also estimated it 
through equation as baryonic disc scale lengths and disc 
rotation velocities for the resulting simulated galaxies were 
given. The test showed an unbiased and tight one-to-one cor- 
relation, with very small dispersion leading to errors < 30%. 
Application to a large volume limited sample of over 11,000 
galaxies from the SDSS yielded a log-normal distribution 
of A for the total sample, interestingly in consistence with 
results from cosmological n-body simulations, presented in 
Hernandez et al. (2007). The above results were then repro- 
duced by Berta et al. (2008), using a larger sample of 50,000 
galaxies also from the SDSS, and a refined version of the A 
estimate given here, for spiral galaxies. 



3 THE SDSS SAMPLE 

The sample of galaxy pairs used in this work comes form 
a study by Park, Gott & Choi (2008, hereafter PGC), us- 
ing data from the SDSS. It is a volume-limited sample of 
galaxies with absolute magnitude M r < —19.5 + 0.51og7i in 
the redshift interval 0.001 < z < 0.5. Since most theoretical 
studies concerning spin distributions present their results at 
z = 0, we limited the sample to low redshifts. 

The nearest neighbour for a given galaxy is found re- 
quiring the following conditions; (1) the neighbour galaxy 
can not be fainter than the target galaxy by more than 
AM r , (2) it must have the smallest projected separation 



across the line of sight from the target galaxy and (3), it 
must present a radial velocity difference less than Vmax ■ We 
choose AM T = 0.5 to include only influential neighbours and 
not merely satellite galaxies. In PGC it was shown that the 
selection of target galaxies having neighbors fainter by more 
than 0.5 mag produces similar results but drastically reduces 
the number of target galaxies as their absolute magnitude 
cut becomes brighter, which reduces the statistical signifi- 
cance of the results. 

To determine Vmax, PGC searched for all neighbours 
with a velocity difference of < lOOOfcms -1 with respect to 
each target galaxy and with a magnitude not fainter by more 
than AM r - When looking at the rms velocity difference of 
the neighbours as a function of the projected separation, this 
remains constant out to 50k" 1 kpc, at 255 and 169 fcms -1 
for early and late type target galaxies respectively. In this 
way, we adopt Vmax = 600 and 400 fcms -1 for early and late 
type galaxies, limits that correspond to about 2.3 times the 
rms values. 

Given that the limiting magnitude of the sample is 
M r = —19.0 + 51og/i, we study only those target galaxies 
brighter than M r = —19.5 + 51og/i, in order to avoid los- 
ing neighbours, we fixed our attention on pairs of relatively 
bright galaxies. For more details of the sample see Park, 
et al. (2008). For the same sample, Choi et al. (2007) have 
determined the exponential disc scale, absolute magnitude, 
velocity dispersion, de Vacouleurs radius and eccentricity for 
each galaxy, assuming a ACDM universe with Qm = 0.27, 
n A = 0.73 and h = 0.71. 

In order to discriminate between elliptical and disc 
galaxies, we used the prescription of Park & Choi (2005) 
in which early (ellipticals and lenticulars) and late (spirals) 
types are segregated in a u — r colour versus g — i colour 
gradient space and in the concentration index space. They 
tested extensively the selection criteria through direct com- 
parison of visually assigned types for a large sample of sev- 
eral thousand galaxies. The specific selection criteria can be 
found in Park & Choi (2005), but essentially select as early 
types, galaxies with red colours, minimal colour gradients 
and high concentration indices. Throughout the paper we 
will consider ellipticals and lenticulars as early type galax- 
ies, and spirals as late type ones. 

In the next section we will examine the correlation be- 
tween the angular momenta of galaxies in pairs, taking into 
account angular momentum magnitude through the A spin 
parameter, and direction through position angle. Later we 
will investigate if the angular momentum of a given galaxy is 
modified by the presence of a neighbour. To conduct these 
inquiries, we obtained four different subsamples according 
to the following criteria. 

Our sample A, used to search for a correlation in the A 
value of spiral galaxy pairs, contains 347 pairs of late type 
late - type galaxies, with 0.001 < z < 0.5, limiting mag- 
nitude of Mr — —19.5 + 51og/i, and Rd measured by Choi 
et al. (2007), and rotational velocities infered from the ab- 
solute magnitude introducing a TF relation (Pizagno et al. 
2007), which is all the information needed to calculate A 
from equation [5] To avoid the problem of internal absorp- 
tion in edge-on galaxies (Unterborn & Ryden 2008; Cho & 
Park 2009), and consequently underestimating rotational ve- 
locities, we limit the sample to spiral galaxies having axis 
ratios b/a > 0.6. 
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To study the correlation of the angular momentum di- 
rection between members of a spiral pair, we constructed 
sample B, composed of 218 pairs, with the same redshift 
range and limiting absolute magnitude as sample A. Sup- 
posing the position angle of observed galaxies to be per- 
pendicular to their angular momentum vector, we have an 
indicator of the angular momentum orientation, albeit with 
a degeneracy between parallel and anti parallel spins. Given 
the difficulty of measuring the position angle for face-on late- 
type galaxies, we used pairs of spiral galaxies with b/a < 0.7. 
This is the reason why we could not use the same sample 
A to search for a correlation in the direction of the angular 
momentum. 

In both samples, A and B, we required that the near- 
est neighbour of a given neighbour galaxy in a galaxy pair, 
should be the target galaxy itself, in this way we search 
for spin correlations for pairs involving only two galaxies, 
the simplest scenario. To test the response of the spin to 
the presence of a companion, we relaxed both this condition 
and the condition of the neighbour galaxy being a late type 
galaxy, which increases the number of galaxy pairs to 3624 
for sample C, which will be used to measure the response 
of the value of A to the distance of the nearest neighbour, 
and to 2037 for sample D with the position angle well deter- 
mined for every galaxy involved, to determine the influence 
of a neighbour on alignment. 

When studying the influence of a nearby galaxy on the 
value of the spin, the distance between galaxies is normalized 
to the virial radius of the neighbour. We define the virial 
radius of a galaxy as the projected radius where the mean 
mass density within the sphere of radius r v t r is 200 times 
the critical density or 740 times the mean density of the 
universe. 



with the relative mass-to-light ratios for early type 
galaxies (ellipticals and lenticulars) twice the same ratio for 
the late types (spirals), y (early) = 2^(late), following Choi 
et al. (2007), who report that the central velocity dispersion 
of early type galaxies brighter than M r — —19.5, is about 
y/2 times that of late types. Since we adopt Qm = 0.27, 
we have 200p c = 2007j/^a/ = 740p, which is almost equal 
to the virialized density p vir — 187r 2 /r2 m (_ffoto) 2 p = 766p, 
in the case of a ACDM universe (Gott & Rees 1975). Fi- 
nally, we introduced the mean value for the density of the 
Universe, p = (0.0223 ± 0.0005)(7L)-2o(/i _1 Mpc)~ 3 , where 
(7L)_2o is the mass of a late type galaxy with M r = —20 
(Park et al. 2008). In this way, the virial radii of galaxies 
with Mr = -19.5, -20.0 and -20.5 are 260, 300 and 350/i _1 
kpc for early types and 210, 240 and 280h~ 1 kpc for late 
types, respectively. 

Having limited ourselves to such relatively massive sys- 
tems, which exclude any dwarf galaxies, also gives confidence 
on the validity of the constant baryon fraction hypothesis 
used in the estimates of A. The above, as small systems 
where substantial mass loss due to the feedback effects of 
star formation probably applies, are not included. 



4 RESULTS 

The stringent selection criteria described above ensure that 
our samples contain only galaxies where our estimates of the 
spin are most reliable, although the samples are reduced to 
a small fraction of the total SDSS field, they remain statis- 
tically significant. 

The results of our A estimates and position angle stud- 
ies appear in this section. Section 4.1 tests the angular mo- 
mentum acquisition mechanism through an exploration of 
the correlations in magnitude and orientation of the angu- 
lar momentum for close pairs of spiral galaxies. In section 4.2 
we investigate the effect on spin magnitude and orientation, 
of the interaction of a spiral galaxy with another galaxy, as 
a function of the separation between them, regardless of the 
morphology of the companion. 

4.1 Spin correlations 

It is reasonable to think that, given the proximity of the 
galaxies in each pair, both galaxies are immersed in the 
same tidal field. The amount of angular momentum gained 
by each galaxy should be similar, proportional to the tidal 
torques exerted by the surrounding material and inversely 
proportional to the mass of each galaxy. This assumes that 
the two components of the present pair have been in close 
association, at least since the phase of angular momentum 
acquisition. We explore this hypothesis with the selected 
sample of pairs of spiral galaxies, calculating spin parame- 
ters using equation [5] and the total mass from the baryonic 
Tully-Fisher relation introduced to obtain equation [5] 

If both member galaxies of a pair were exposed to the 
influence of the same tidal field, their product AAf should 
show a clear and tight 1 to 1 correlation. In Fig. 1 top panel 
is shown the spin mass product of nearest neighbour galaxies 
M2A2, as a function of the same product for target galax- 
ies M1A1, in units of solar masses for the 347 pairs of late 
type - late type galaxies extracted from our sample A. In 
the scenario where both galaxies gained their angular mo- 
mentum through interactions with a constant tidal field, we 
should see a clear correlation, but the correlation between 
the products for neighbouring galaxies in this case is low, 
with a correlation index of r 2 — 0.247. 

To quantify the low level correlation seen in the top 
panel of figure 1, we produced a large number of equiva- 
lent mock catalogues of spiral-spiral galaxy pairs through 
Monte Carlo simulations. Each, with a total of 347 pairs, 
was obtained through sampling a fixed distribution function 
for both of the numbers which describe the members of an 
observed pair. A random value of M and A is picked for each 
member of a sampled pair, until one obtains as many mock 
pairs as present in our sample A. 

Given the restriction in the absolute magnitude im- 
posed when selecting the observed pairs, the first condition 
to find the nearest neighbour, a constraint is imposed upon 
the difference in mass within the members of a given pair. 
The mass input for the test was chosen observing this restric- 
tion, using a constant density distribution with upper and 
lower limits for the mass as imposed on the SDSS pairs sam- 
ple. The values of A were taken from the log-normal distri- 
bution obtained empirically in Hernandez et al. (2007) with 
parameters Ao = 0.0394 and a\ = 0.509, for a larger sample 
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Table 1. Correlation indexes for Sample A galaxies and mock 
catalogues, showing for the latter (r 2 ) ± a values. 
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Figure 1. AiMi of target galaxies against A2M2 of nearest neigh- 
bour galaxies, in units of solar masses, for 347 pairs. Top panel: 
Pairs of sample A with a r 2 = 0.247. Bottom paneLfor compari- 
son, result from a Monte Carlo sample with no inherent correla- 
tion beyond what is imposed by the magnitude selection criteria 
in the observed sample, with r 2 = 0.045. 



of SDSS galaxies. A mock sample is shown in Fig. 1 bottom 
panel, for comparison with the sample A of galaxy pairs from 
the SDSS, showing a correlation index of r 2 = 0.135. In this 
particular case, the correlation index is lower for the mock 
sample, perhaps the result of an statistical fluctuation. To 
quantify this effect we measured the correlation index for 
one thousand mock catalogues of pairs, all with 347 pairs, 
and found a mean value of (r 2 ) = 0.129 and, given that the 
collection of synthetic correlation indexes is well described 





r 2 


Sample A 


mock catalogues 


M2A2 


vs AfiAi 


0.247 


0.129 ±0.063 


A 2 


vs Ai 


0.202 


0.0056 ± 0.0604 



by a normal distribution with a spread of 0.063, it places 
the correlation between the SDSS pairs about 1.9a above 
the expected for null intrinsic correlation between A1M1 and 
A2M2. Even though the correlation for the sample of SDSS 
galaxy pairs seems weak and the correlation index is low, the 
correlation is significantly stronger than that obtained for 
null intrinsic correlation mock samples. We must take into 
account that the correlation extends over 1.5 magnitudes 
in the MX product, and that this relation is dominated by 
galaxy pairs with large separations. 

Due to the restriction in magnitude imposed when 
defining the sample of observed pairs, the difference in mass 
between each member of a pair is small, and allows us to 
search for a direct correlation between Ai and A2. The result, 
using the same sample A, shows no difference with the pre- 
vious test, with a r 2 = 0.202, a weaker value but still larger 
than what is obtained with the results of a Monte Carlo test 
with no intrinsic correlation at all; with (r 2 ) = 0.0056 and 
a — 0.0604. It appears likely that the correlations shown 
in Table 1 have their origin in a common tidal field having 
being at least partly responsible for imprinting the angular 
momentum to both members of the pairs studied. However, 
the fact that this correlations are weak, points to other an- 
gular momentum acquisition and modification mechanisms 
having being at work, e.g. mergers. 

So far we have centred our attention on the magnitude 
of the spin parameter, but we have not yet examined the 
direction of the angular momenta in the pairs. 

For the analysis of the angular momentum direction, 
we use sample B. We measured the difference in the posi- 
tion angle for the galaxies of each pair, restricted to vary in 
the range < < 90, given our impossibility to discrimi- 
nate between parallel and anti-parallel spins. Once obtaining 
for each pair, following the method proposed by Yang et 
al. (2006), we count the total number of pairs, N (0), for a 
number of bins in 0. Next, we construct 100 random sam- 
ples in which we randomized the orientation of the galaxies, 
respecting the selection criteria of the real sample, and we 
computed the average number of pairs {Nn (0)), as a func- 
tion of 0. By construction, the random samples have the 
same selection effects as the real sample of pairs, so any 
significant difference between N (0) and Nr (0) reflects a 
genuine alignment between galaxy pairs. 

To quantify the strength of any possible alignment, we 
define the normalized pair count as 



fpairs (0) 



N(0) 



(7) 



(JVh(0)>- 

For a complete absence of alignment we should obtain 
fpairs — 1. To assess the significance of the deviation of the 
normalized pair count form unity, we use o~r (0) / (Nr (0)}, 
where or (0) is the standard deviation of Nr obtained form 
the 100 random samples. 

The result of applying this technique to our sample B is 
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Figure 2. Normalized probability distribution of the difference 
in the position angle, ©, for the pairs in sample B, dotted line 
representing f pa ira = 1] and the error bars give the value of 

°r (e)/(N R (&)). 

presented in Fig. 2, where we show f pa irs for the 255 pairs in 
the sample, and the value of or (6) / {Nr (©)) as error bars. 
No intrinsic alignment is evident, as there are no significant 
deviations from f pairs = 1, represented by the dotted line, 
which is consistent with the result obtained with the A mag- 
nitudes of the pairs; a very weak correlation between the 
spins of neighbouring spiral galaxies, consistent with pre- 
vious studies finding little correlation between the position 
angles of spiral-spiral pairs (Wang et al. 2009), which hence 
appear to have formed as a pair mostly after the phase of 
angular momentum acquisition was over. In this particular 
case, the sample is dominated by pairs with large separa- 
tions; i. e. d > R V ir, as we will see in the next subsection. 

The low correlation between the spin of neighbouring 
galaxies can be explained by the early acquisition of angular 
momentum through primordial torques exerted by the sur- 
rounding tidal field on the pair (Barnes & Efstathiou 1987; 
Navarro, Abadi & Steinmetz 2004) , which induces an initial 
inherent correlation. In this sense, the correlations found 
(reported in table 1) constitute an important direct empiri- 
cal confirmation of the validity and importance of primordial 
tidal torquing in contributing to galactic angular momentum 
acquisition. The effects of different subsequent mechanisms 
such as strong interactions with neighbour galaxies (Hern- 
quist 1993), delayed mass aggregations events (Kitzbichler & 
Saurer 2003; Pierani et al. 2004) and minor mergers (Gard- 
ner 2001; D'Onghia & Navarro 2006) are to then blur the 
initial correlations, to the point where no evident alignment 
correlation remains, highlighting the effect of later processes. 
In the next subsection we will focus on the effect of interac- 
tions in the evolution of the spin. 

4.2 Interacting galaxies 

In this subsection we investigate if there is a response to 
the interaction between neighbouring galaxies in the A spin 



parameter, as a function of the separation distance between 
the members of the observed pairs. In this case we are only 
concerned with the value of A for late type target galaxies, 
irrespective of the A characteristics of the nearest neighbour. 
This allowed us to increment the number of systems studied, 
through the inclusion of spiral target galaxies having early 
type closest neighbours, for the third of our tests. 

Recently Park & Choi (2009) reported a dependence of 
the absolute magnitude of the target galaxy on the nearest 
neighbour separation, where the galaxy luminosity in the red 
decreases as the separation distance decrease. Given that in 
the calculation of A we use M r to assign a Vd value for equa- 
tion [5j we limited the sample to galaxies having magnitudes 
in the range —19.5 < M r < —20.5, to focus our study on 
the response of A not merely as a consequence of the depen- 
dence found by Park & Choi (2009) on magnitude, but as a 
directly consequence of the influence of a neighbour galaxy 
on the value of A. Taking only a narrow range in M r limits 
the effect on the estimates of A of a correlation between M r 
and the nearest neighbour distance, to better isolate and 
asses the effects of the distance to the nearest neighbour, on 
the spin of the target galaxy. 

The top panel of figure 3 shows A values of the 3624 
target late type galaxies of sample C, as a function of the 
distance to their closest neighbours, normalized to the virial 
radius of the neighbour galaxy, calculated using equation 
[U The sample is divided into 7 bins, where the median A 
values are shown with their dispersion presented as thin 
error bars and the uncertainty represented by thick error 
bars; this convention will be followed for the next figures. 
For separations larger than the virial radius of the neigh- 
bour galaxy, the mean value of A appears constant, but as 
soon as the distance becomes smaller than R V i r , the me- 
dian value of A starts to decrease, a clear indication of the 
interaction. Galaxy harassment and various angular momen- 
tum loss mechanisms, such as dynamical friction (Hernquist 
1993), might begin to operate during the interaction, even 
at large distances, as the galaxies first cross into their virial 
radii (PGC; Park & Choi 2009). 

In going to the first bin, a substantial increase in the dis- 
persion of the inferred values of A is evident. However, for 
such small separations, systems within 0.02 of their virial 
radius, a strong interaction is ongoing. This implies the sys- 
tems are heavily disturbed and strongly out of equilibrium; 
see e.g. the numerical simulations of galaxy collisions of Her- 
nandez & Lee (2004), where fluctuations in potential and ki- 
netic energies of the total system are assessed as a function 
of the pair separation. This invalidates the assumptions go- 
ing into the simple formula used to estimate A, making the 
values of this quantity in the first bin of figure 3, top panel, 
useful only as indications of strong interactions at small dis- 
tances, and not as indicative of total halo A. 

The complementary study using the difference in the 
position angle O, to account for the relative orientation of 
the galaxies within a pair, is shown in Fig. 3 bottom panel 
using sample D, where Q is plotted against the separation 
distance between neighbouring galaxies, normalised by the 
virial radius R V i r of the neighbour galaxy. The sample is 
divided into 7 bins, where the median values are shown 
with their dispersion presented as thin black error bars, as 
in the preceding plot. If the orientations of the galaxies be- 
tween neighbours were isotropic, we would expect a median 
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Figure 3. Top panel: A value for the 3624 target late type from 
sample C as a function of the separation distance to their near- 
est neighbour, normalized by the virial radius of the neighbour 
galaxy, R v i r bottom panel: Difference in the position angle © for 
the 2037 late type pairs from sample D as a function of their 
normalized separation distance d/R v i r 

O value of 45, showed as a broken line in the plot. As can 
be seen, the median value of O decreases as the separation 
distance decreases below the virial radius of the neighbour 
galaxy, meaning that the galaxies tend to be aligned (or 
anti-aligned) as their neighbours get closer to them. 

Summarising, we have shown that the angular momen- 
tum of spiral galaxies, not only in its orientation but also in 
its magnitude, is affected by the presence of a companion, 
once the galaxies are within the virial radius of their neigh- 
bour. Once galaxies cross into their virial radii, a slight but 
measurable decrease in A begins to operate. A corresponding 
increase in alignment also appears, but only in going to very 
close systems where the out-of-equilibrium situation makes 
determining values of A unfeasible. 

In a previous work (Cervantes-Sodi & Hernandez 2009) 
using a reduced sample of very well studied spiral galax- 



ies, two of us showed how the star formation rate has a 
marked dependence on A, it decreases when the value of A 
increases. If such tendency remains in perturbed systems, 
it could explain the enhanced star formation present in in- 
teracting galaxies, as due to the decrease of the spin. Some 
studies (Lambas, Tissera, Alonso & Coldwell 2003; Ellison, 
Patton, Simard & McConnachie 2008) based on observa- 
tional samples coming from large galaxy surveys, even show 
a gradual increase of the star formation as the distance be- 
tween the interacting galaxies decreases, result completely 
compatible with our finding of the decrease of the value of A 
in such circumstances. More recently, Park & Choi (2009), 
studying the effects of galaxy interactions on galaxy proper- 
ties, reported the enhancement of the star formation activity 
for target galaxies within the virial radius of the neighbour 
galaxy, when this neighbour is a late type galaxy. 

We explore the effects of the galaxy type of the nearest 
neighbour by splitting sample C in two according to the 
neighbours morphology, and plotting the value of A as a 
function of the normalised distance between neighbouring 
galaxies d/R V i r , in Fig. 4 left panels. The top panel shows 
the behaviour when the nearest neighbour is a late type 
galaxy (2073 galaxies) and bottom panel when it is an early 
type one (1551 galaxies). The complementary plots showing 
the difference in the position angle splitting sample D in two 
according to the neighbours morphology, are shown in the 
right panels of Fig. 4; top panel for nearest neighbours being 
late type (1325 galaxies) and bottom being early type (712 
galaxies). If we look at the case when the neighbour is a 
late type galaxy, we see that within the errors, no effects are 
apparent when the normalised distance is larger than 1, with 
a slight decrease in the values of A as the normalized distance 
drops below lR V i r and with an increase in the alignment or 
anti-alignment between their position angles. The case when 
the neighbour galaxy is an early type one is more erratic and 
the trend less evident, in both the magnitude of A and the 
spin alignment. It is interesting to note that for the first 
bin, where the interaction is ongoing, we see much more 
significant disturbances in the value of A, when both galaxies 
are spiral, probably due to the strong hydrodynamical effects 
of the interaction of two gaseous disks, which are absent 
when one of the two systems is of early type. 

As an extension to the work presented here, it would be 
desirable to probe particularly high density environments, 
to test how the situation tends to the picture seen at cluster 
scales (see Park & Hwang 2009 for environmental effects on 
cluster galaxies). Within the hierarchical picture of structure 
formation, evolutionary timescales decrease with increasing 
ambient densities, which in turn also increase with time. In 
this way, one could hope to recover an evolutionary sequence 
for any given process, from a comparative study at zero red- 
shift, extending over a large ambient density range. Also, a 
fuller 3D alignment study might yield more extensive infor- 
mation. One must be cautious, e.g., on cluster scales, results 
of 2D alignment studies are sometimes modified when more 
extensive 3D analysis are performed. 



5 CONCLUSIONS 

Using pairs of spiral galaxies of comparable size, we searched 
for correlations between their spins. Using the spin param- 
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Figure 4. Left panels: A value of target late type galaxies from sample C as a function of the normalized separation distance to their 
nearest neighbour, d/R v i r ; top panel: nearest neighbour being a late type galaxy, bottom panel: nearest neighbour being an early type 
galaxy. 

Right panels: difference in the position angle for galaxies from sample D as a function of the normalized separation distance to the 
nearest neighbour; top panel: nearest neighbour being a late type galaxy, bottom panel: nearest neighbour being an early type galaxy. 



eter A to account for the magnitude of the angular momen- 
tum, and the position angle for its direction, we found a 
weak but statistically relevant correlation between the spin 
magnitude of neighbouring galaxies, but where unable to 
detect any tendency for alignment. 

If we adopt the simplest version of the tidal torque the- 
ory to explain the acquisition of angular momentum, we 
would expect a strong correlation for the spins of spiral 
galaxies in a pair, if they had formed under the influence 
of the same tidal field. Our results imply the dominant pres- 
ence of two complementary effects, the late formation of the 
galaxy pair, after the epoch of A acquisition was completed, 
and the blurring of initial conditions due to the clumpy and 
irregular mass accretion of relatively extended minor merg- 
ing events. 

In going to the trends exhibited by the sample as a 



function of galaxy separation, we can clearly see the effects 
of galaxy-galaxy interaction. These appear as soon as the 
galaxies cross into their virial radii, leading to a gradual 
decrease in the values of A. 



Our results do not diminishes the importance of the 
torquing at early stages of galaxy formation, but give us 
an insight into the important role of later interactions in 
the overall evolution of the total angular momentum of the 
galaxies, as shown by the results obtained from interacting 
systems, where the spin is visibly perturbed by the presence 
of a nearby companion. 
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